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We present a mean-field theory for magnetic field driven transitions in dipolar coupled gadolinium 
titanate Gd2Ti2C>7 pyrochlore system. Low temperature neutron scattering yields a phase that can 
be regarded as a 8 sublattice antiferromagnet, in which long-ranged ordered moments and fluctuating 
moments coexist. Our theory gives parameter regions where such a phase is realized, and predicts 
several other phases, with transitions amongst them driven by magnetic field as well as temperature. 
We find several instances of local disorder parameters describing the transitions. 



I. INTRODUCTION 

The pyrochlore S=7/2 Heisenberg spin system 
Gd2Ti207 is currently very popular. In addition to mag- 
netic ordering at 1 Ki which is not expected for a Heisen- 
berg pyrochlore system^ it displays magnetic field-driven 
phase transitions that are intriguing, and have been as- 
cribed to the competition between dipole-dipole (d-d) 
interactions and the Zeeman energy on frustrated and 
undistorted cubic systems* 3 " The exchange energy has 
been argued to be relatively small, due to the compact 
f-shell of Gd, but nevertheless plays an important role 
in lifting the degeneracy between various possibilities, as 
we show in this work. 

Earlier examples of similarly low energy scales are 
known, but on the diamond lattice^ and the gadolinium 
gallium garnet^ thus lacking the complexity of the py- 
rochlorcs. Theoretical interest in d-d coupled systems is 
quite old, 6 focussing on its long ranged nature. 

Gd2Ti2C>7 is intriguing, first by the magnetic order it 
displays at low temperatures and zero magnetic field, as 
observed by neutron scattering^ The magnetic order has 
been assigned a propagation vector q = {tt, tt, it} = tt at 
50 mK, and the interstitial moments^ are zero on average, 
thus showing a partially-ordered state despite the very 
low temperature^ This structure is consistent with the 
correlations observed at higher temperatures in ESR. 9 
This turns out to be in direct conflict with previous theo- 
retical works that assumed a nearest neighbor Heisenberg 
system with d-d interactions. The expansion of the free- 
energy to quadratic order in the order parameter near 
the ordering temperature, T/v, tells us that all the states 
with a propagation vector along the (111) direction be- 
come unstable simultaneously.' 1 ' Such a degeneracy is in 
fact not quite exact; refined numerics and analytical work 
show that the q = tt state is very slightly preferred^ 
Nevertheless, below T/v the fourth order term in the free 
energy expansion becomes important and favors a q = 
state^i which has been confirmed by a real-space mean- 
field theory employing a four sublattice decomposition 
of the pyrochlore latticed From these works, however, it 
is not clear as to how this q = tt state can be realis- 
tically and robustly realized. We find that this requires 
going beyond the nearest neighbor exchange model stud- 



ied previously. 

The second particularly interesting feature of gadolin- 
ium titanate is the presence of magnetic-field-driven 
phase transitions^ For cubic systems one expects exotic 
phases connected by continuous transitions, as compared 
to well studied first order spin-flop type transition on uni- 
axially distorted systems, such as MnF2. Although the 
four sublattice mean-field theory is in conflict with the 
neutron scattering results, it gives such field-driven tran- 
sitions with critical fields in rough agreement with those 
observed on a powder sampled 

The aim of the present article is to show that a 
model that includes exchange energies beyond the near- 
est neighbors not only explains the magnetic structure 
observed by neutron scattering at low temperatures but 
also exhibits magnetic-field-driven transitions whose fea- 
tures can be checked by further experimentation. We 
present the results of the expansion of the free-energy 
near T/v together with a low temperature mean-field the- 
ory using "hard-spins" (that describes the fixed length 
constraint well) that is, hence, valid at all temperatures, 
although limited by the sublattice structure imposed at 
the outset. The earlier analysis 3 imposed a 4 sublattice 
order, here we go to the next level of description, namely 
a 8 sublattice order, in order to accomodate the q = tt 
state. This allows us to find other stable phases. 

We begin by summarizing the results of the "hard- 
spin" mean-field theory. The 8-sublattice system con- 
sists of two tetrahedra (see Fig. [TJ, the second one be- 
ing obtained by a translation of the first one along the 
primitive vector (110). There are 4 other such equivalent 
directions that would give equivalent results, in particu- 
lar with regards to the propagation vector tt. Without 
loss of generality, we consider only one of them in the 
following, namely (111). 

II. ZERO MAGNETIC FIELD PHASES 

The pyrochlore lattice sustains various phases. We 
have alluded to the q = phase previously, that we 
denote A in the present work. It has a 6-fold degener- 
acy according to the three planes of the cubic structure 
and the time-reversal symmetryi^ii In addition there are 
three other phases that we call B, C and D, that break 
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translation invariance. 

At finite temperatures, and for different (exchange) pa- 
rameters, we find a distinct phase B, that can be de- 
scribed as a q = 7r state. In this phase, the interstitial 
sites have, on average, no magnetic moment. The spins 
of the Kagome planes, on the other hand, are ordered in 
a 120 degree structure; each of them being parallel to the 
opposite edge in order to minimize the d-d interactions. 
This is the phase that has been found experimentally in 
Gd 2 Ti 2 7 . 

There is also a phase that we call C which does not 
have a simple description: the interstitial spins (1 — 1') 
are parallel to each other and the other spins break trans- 
lational symmetry and are not coplanar. There are 12 de- 
generate such states. It seems plausible that the C-type 
phases are in reality the 'projection' onto the 8-sublatticc 
of incommensurate states. 

In addition to these phases which all have a finite de- 
generacy, we have found a phase, D, which is particularly 
interesting in that it has a continuous degeneracy. We 
emphasize here that we actually start from a strongly 
frustrated system, whose huge degeneracy is lifted by 
the d-d interactions. It is therefore surprising to ob- 
tain a continuously degenerate ground state in a wide 
range of parameters when further neighbor interactions 
are switched on. In the D phase, the spins of the Kagome 
planes are identical to those of the B phase. The inter- 
stitial moments, however, are finite. They may point in 
any direction in a plane parallel to the Kagome planes, 
giving a large degeneracy, and they are antiparallel to 
each other from plane to plane. This is actually the zero 
temperature analogue of the B phase, but occurs at low 
temperatures. 
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FIG. 1: On the left, the eight sublattices used in the mean- 
field calculation. The A state (q=0) is shown ( six fold de- 
generate). On the right, the spins of the B state (degeneracy 
2) belong to the Kagome planes; the sublattices 1 and 1' have 
no magnetic moment in average. The D state is identical to B 
except that the magnetizations of the sublattices 1 and 1' are 
finite and opposite from one Kagome plane to the next (con- 
tinuous degeneracy). C does not have any simple description 
(see text). 



To find those phases, we have considered the Hamilto- 
nian defined on the pyrochlore lattice 

H = — ^ JjjSj ■ Sj — g[iB ^ S; • H 

ij 
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where Si is a quantum spin operator for spins S = 7/2 on 
site i. Jij is the Heisenberg exchange between the neigh- 
bors: we are in fact considering the first neighbors at a 
distance a = 3. 59 A (J), the second (J 2 , distance \/3ao) 
and third neighbors (J3, distance 2ao). We note that 
there are in fact two types of third neighbor exchanges 
according to the crystal structure. The difference of them 
being another small parameter, we will neglect it in the 
following. Let us first describe the "k-space" mean-field 
theory. We have studied the stability of the paramag- 
netic state towards a state with a modulation vector k, 
extending RefQ to include J 2 and J3. By expanding 
the free-energy near Tjy, we get a 12 x 12 matrix, its 
lowest eigenvalue determines T/v and the corresponding 
eigenvector describes the ^-modulation of the state. We 
find that the degeneracy along (lll)iil for J 2 = J 3 = 
is actually weakly lifted when the dipolar sums include 
a larger number of neighbors compared with previous 
works, as previously noticed^ The selected mode is at, 
or very close to, q = it (see Fig. New numerical work 
using the Ewald summation technique has indeed clearly 
shown that this mode is selected^ Nevertheless in order 
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FIG. 2: Quasi-degeneracy of the lowest eigenvalue as function 
of q along (111) for J2 = J3 = (see the horizontal scale). 
When L (the number of neighbors in the (100) direction in- 
cluded in the dipolar sums) is increased, the curves converge 
to a smooth curve which minimum is at, or very close to 
q = 7T. 

to robustly lift what remains a quasi-degeneracy, it is im- 
portant to include J 2 and J3. As soon as J 2 < 0, q = n is 
robustly selected. The corresponding eigenvector tells us 
that it corresponds to the B state. There are two other 
regions of the phase diagram where the A phase and an 
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incommensurate phase are preferred. The phase diagram 
of the first instability is given in dashed lines in Fig. [3] 
Since the above approach can only give the first insta- 
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FIG. 3: Phase diagrams showing the first unstable modes of 
the paramagnetic phase at Tat (dashed line, small letters, " i 
" for incommensurate) and from the 8-sublattice mean-field 
theory (solid line, capital letters). In the inset, the phase 
diagram as a function of temperature is shown for Ji — 0.1J. 
Tn is the temperature to enter the B phase at intermediate 
J 3 . J = 0.4 K. 

bility that is encountered on cooling, i.e. at Tjy, nothing 
is known about the lower temperature behavior of the 
system. We therefore proceed to a real-space mean-field 
theory following Ref. y, but enlarging the unit cell to 8 
sublattices, as already specified. 

Within the mean-field theory, the magnetizations of 
the sublattices are given by (S a ) = Sh a Bs(f3h a S) with 
the definitions of the local fields, ft," = gfiBH a — 
Jab \&b) ' w here &s is the Brillouin function, and 
(3 is the inverse temperature. J^f couples the sublattices 
a and b, running from 1 to 8. A straightforward way 
of solving the problem consists of iterating numerically 
these self-consistent equations, starting from many ran- 
dom configurations and selecting the set of states with 
the lowest free-energy. This leads to the phase diagram 
given in Fig. [21 (solid lines). The phases A, C and D 
described above are found at zero temperature in quite 
a narrow and physically reachable region of parameters. 
However, there is no sign at zero temperature of the B 
phase found in the mode analysis. 

In Fig. |3 we give the labels for the sequence of phases 
on increasing T prior to reaching the paramagnetic phase. 
Details of the transition temperatures for the 8-sublattice 
problem are given in the inset of Fig. [3] for a given J2 = 
0.1J (solid lines). For J2 = J3 = 0, as another example, 
the A phase is definitely chosen below Tn, in agreement 
with RefslMrj. The transition to the B state at Tjy, as 
suggested in the mode analysis given above (fig. EJ, is 
in fact immediately followed by a transition to the A 
phase when the temperature is decreased (so that on the 
scale of fig. 21 the B phase is indistinguishable at the 
PM-A transition). To robustly select the B phase in a 



wide range of temperatures (as experimentally observed) 
it is necessary to include other interactions, the simplest 
assumption being the next nearest Heisenberg couplings. 
Indeed when J3 is increased or J2 decreased, the region 
where we have the A-B succession of phases widens (see 
inset of Fig. |3J, before the onset of the other phases. For 
these other phases, the C-B transition line appears to be 
first-order, whereas D-B is second-order. 

We can now compare in more details these results with 
those of the fc-space mean-field theory, i.e. the stability 
analysis of the paramagnetic state. For large portions of 
the phase diagram (where the last capital letter coincides 
with the small letter), the two approaches give the same 
result. There is, however, a region of parameters (e.g. 
large J3) where the first unstable mode is incommensu- 
rate, while for the same parameters wc find a transition 
from the paramagnetic state to the B state in the 8- 
sublattice calculation. The 8-sublattice calculation can 
not capture the first transition to the incommensurate 
state because of the sublattice decomposition. On the 
other hand, such a calculation respects the fixed length 
constraint of the spins at low temperatures. Within this 
point of view, two scenarios may take place when the 
temperature is lowered. In both of them there is first a 
transition from the paramagnetic state to the incommen- 
surate state. Then the incommensurate state may exist 
down to zero temperature and the occurrence of the B 
state is a pure artifact of the sublattice calculation. Al- 
ternatively, the B state is stabilized at low temperature. 
In this case there must be a phase transition between the 
incommensurate state and the B state when the temper- 
ature is lowered. In order to rule out one scenario, one 
would need to study the stability of the phases at low 
temperature bypassing the limitations of the sublattice 
decomposition. 



III. MAGNETIC-FIELD-DRIVEN 
TRANSITIONS 

A magnetic field is an interesting probe of these phases. 
We highlight here the main results on all four phases since 
we believe that the A, C and D phases might be useful for 
other materials. The A phase gives rise to multiple phase 
transitions, as previously reported for T=0£ We note, 
however, several new elements arising from the effect of 
finite temperatures. In Fig. we give the example of 
the complete phase diagram when the field is along (110). 
For small enough fields, a raise in temperature drives 
the system through two phase transitions, a result which 
also holds for the fields along the other crystallographic 
directions, (100) and (111), with, in the latter case, a 
first transition which is weakly first-order. 

Although the two successive phases A^ for a field along 
(110) are separated by a transition line (see Fig. QJ, they 
can not be distinguished by a different geometrical broken 
symmetry. They basically bear the same ad broken sym- 
metry, but none of the other geometrical symmetries is 
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FIG. 4: Phase diagram (H,T) for a magnetic field parallel 
to (110), showing the broken symmetries (in brackets). A n 
represents the A phase under field. The B phase occurs in a 
very narrow region near Tjv at H = (see text). 



broken in either phase. The effect of thermal fluctuations 
turns out to be crucial in distinguishing them. Indeed, 
at any finite temperatures, this transition is associated 
with the occurence of a disorder parameter - a quantity 
which is zero in the ordered phase and non-zero in the 
disordered phase, which we now describe. When we look 
at the value of the magnetic moment (smaller than 5* as 
soon as T > 0), it appears that in the ordered phase 
all the magnetic moments of the four sublattices are the 
same, while they are different above the first transition 
(Fig. [5). We can thus define a disorder parameter by 
the difference of two of these moments (inset of Fig. [5J . 
Moreover, we have found that the disorder parameter has 
a critical exponent ~ 4/3. This is in contrast with the 
other linear order parameters which follow the usual 1/2 
exponent of the Landau theory. One remarkable feature 
to note is that our disorder parameters are local. This 
is quite unlike familiar disorder parameters that arise in 
theories with duality that lead naturaly to highly nonlo- 
cal disorder parameters, such as are known for the 2-d 
Ising model. 

We note that these transitions are quite different from 
the classical spin-flop transitions of the uniaxial antifcrro- 
magnets. In addition to those unconventional mean-field 
exponents, the response of the spins themselves is worth 
mentioning. When a field along (100) is increased, for 
instance, the spins 1 and 2 of the A phase (or 3 and 4) 
which are perpendicular to each other at zero field, re- 
main exactly perpendicular while the sum of both spins 
twist towards the field. This is not obvious and in par- 
ticular can not be understood by considering partial cou- 
ples, since it is a real 4-sublattice-coupled system. 

The effect of a magnetic field on the B-phase is to in- 
duce back a net magnetic moment on the previously zero 
moment interstitial sites. At higher fields along the (111) 
or (100) directions, there is a unique transition to the 
paramagnetic phase with a merging of the two degenerate 
states (Fig. [5J). For the (110) direction, however, there is 
a reentrance of a less symmetric phase when the magnetic 
field is increased. This phase is 4-fold degenerate and 
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FIG. 5: Local magnetic moments averaged by the thermal 
fluctuations as function of temperature. Above the first tran- 
sition, we see a difference in the moments of the two sublat- 
tices that are inequivalent because of the field. The difference 
is a disorder parameter for one transition, with a critical ex- 
ponent numerically close to 4/3. 

breaks the mirror plane symmetry. At larger fields the 
2-fold degenerate state is recovered before reaching the 
transition to the paramagnetic state. In total, there are 
three main distinct transitions when the field is increased 
from zero because the critical field of the paramagnetic 
transition depends very weakly upon the direction of the 
field. 

Regarding the C or D phases, we have to describe what 
happens to the interstitial moments before the system en- 
ters the B-phase and undergoes the transitions described 
above. For the D-phase, the continuous degeneracy is 
preserved with a field along (111), until the B-phase is 
reached. For the other directions, the degeneracy is lifted 
at an infinitesimal field by a spin-flop transition for the 
interstitial moments. For the C-phase, the situation is 
very similar (see Fig. [SJ). 



IV. COMPARISON WITH EXPERIMENTS 

We now compare our theoretical calculations with ex- 
perimental results on Gd2Ti20y. The success of this ap- 
proach is to predict the existence of the q = n phase in 
a wide range of temperatures and parameters. Although 
it is not the ground state of the problem (since only A, C 
and D phases are obtained at zero temperature), the B 
phase may be stabilised at very low temperatures, as ob- 
served experimentally. This is obtained with parameters 
J r 2 and J3 typically of the order of 0.1J which are phys- 
ically reasonable. Next nearest neighbor Heisenberg in- 
teractions may indeed take place in these materials, espe- 
cially if we think of an exchange mechanism in terms of a 
magnetic exchange between the f-electrons and the more 
extended d electrons which may carry the spin polariza- 
tion at distances larger than the nearest neighbors^ We 
give the phase diagram for typical parameters in Fig. \7\ 
which reproduces the number of field-driven transitions 
observed experimentally. We have two free parameters, 
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FIG. 6: Evolution of the degeneracies as a function of an 
external magnetic field along the main crystallographic di- 
rections. The broken symmetries and the symmetries that 
connect the distinct states are also given (T is the translation 
operation, a a mirror plane, and C3 a three- fold rotation). 
Dark shading indicates continuous degeneracy of the state. 
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FIG. 7: Phase diagram (H,T) for a magnetic field parallel to 
(lj_0), showing the various phases and the broken symmetries 
(in brackets). Djj represents the D phase under field. Insert: 
Spin susceptibility for H | (110), convoluted with a gaus- 
sian function (solid line) and powder-averaged susceptibility 
(dashed line). 



J 2 /J and J3/J, chosen typically as —0.125,0.1, and J is 
then fit to the measured high-temperature susceptibility. 
The values of the critical fields measured at low temper- 
atures (3T, 6T and 7T) are then well reproduced. Spin 
susceptibilities are given in the inset of Fig. [7] Note that 
the powder-averaging susceptibility is almost featureless 
and does not account for the complexity of the transi- 
tions. 

Also interesting are the ESR results in the paramag- 
netic phase ^ In addition to anisotropic properties, the 
ESR signal shows two lines, revealing the existence of 
two inequivalent sites, though since no structural distor- 
tions have been reported, the four sites of the unit-cell 
should be equivalent. We note for example that distor- 



tions that would strengthen the magnetic interactions in 
the Kagome planes would naturally help to select the 
B phase. Even in the absence of distortions, the ap- 
plied magnetic field (which is necessary to perform the 
ESR studies) breaks some symmetries of the lattice. In 
this respect, and although the sublattices are equivalent 
at zero magnetic field, we have found different averaged 
magnetic moments on two different sites under magnetic 
field. Therefore, we naturally expect two lines to appear 
in the high-temperature ESR spectrum. Note in partic- 
ular that the difference of magnetic moments is precisely 
what we found to be a disorder parameter in one tran- 
sition at low temperature (fig. ®. In any case, it is 
an interesting example where ESR would not probe the 
properties of the zero-field state. This is of course ex- 
pected from symmetry considerations, but we have found 
here a simple example of a system where it does occur. 

The phase diagram as a function of the model pa- 
rameters (fig. OU may be relevant to other compounds 
as well. The chemical replacement (of Ti by Sn for in- 
stance), or an external pressure could well drive the sys- 
tem into either of the other phases because of the modi- 
fications in the magnetic exchange paths, and hence the 
relative strengths of J, J 2 , and J3. A different behav- 
ior has indeed been found experimentally in Gd 2 Sn 2 07 
at low temperatures^ Although this system orders at a 
similar temperature, the moments were suggested to be 
perpendicular to the local (111) directions on the basis 
of Mossbauer measurements*^ Such a magnetic order- 
ing, if confirmed by neutron scattering, would be more 
compatible with the A phase. 



V. CONCLUSION 

In summary, we have found several phases and phase 
transitions in the dipolar pyrochlore lattice by taking into 
account the magnetic exchange beyond the nearest neigh- 
bour Heisenberg coupling. One of them is precisely the 
B phase of Gd 2 Ti 2 07 seen in neutron scattering^ We 
find that it is not the ground-state, but a finite tempera- 
ture state. We obtain, as ground states, the A, C or the 
degenerate D phase. When a magnetic field is applied, 
we find that the B-phase undergoes three transitions, as 
appeared in specific heat measurements. We have pre- 
dicted the corresponding magnetic structures and the 
broken symmetry phases that could be checked by neu- 
tron diffraction experiments. Sound velocity and absorp- 
tion studies, as well as calorimetric studies with aligned 
crystals in magnetic fields should shed light on the nature 
of the phases predicted here. 
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